Static dielectric properties of carbon nanotubes from first principles 
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We characterize the response of isolated single- (SWNT) and multi-wall (MWNT) carbon nan- 
otubes and bundles to static electric fields using first-principles calculations and density-functional 
theory. The longitudinal polarizability of SWNTs scales as the inverse square of the band gap, while 
in MWNTs and bundles it is given by the sum of the polarizabilities of the constituent tubes. The 
transverse polarizability of SWNTs is insensitive to band gaps and chiralities and is proportional to 
the square of the effective radius; in MWNTs the outer layers dominate the response. The transverse 
response is intermediate between metallic and insulating, and a simple electrostatic model based 
on a scale-invariance relation captures accurately the first-principles results. Dielectric response of 
non-chiral SWNTs in both directions remains linear up to very high values of applied field. 
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PACS numbers: 73.63.Fg; 77.22. Ej; 85.35.Kt 

Carbon nanotubes attract a lot of scientific interest due 
to their unique and versatile electronic and mechanical 
properties, suitable for a wide range of applications. Nan- 
otubes have different electronic properties, determined in 
the zone-folding scheme by the chiral vector: armchair 
(to, to) nanotubes are ID metals, and zigzag (to, 0) nan- 
otubes are semiconductors, with (almost) vanishing gaps 
for to = 3n. Synthesis and separation of specific nan- 
otubes remains a central challenge. Variations in chiral- 
ity and size influence dielectric properties, which in turn 
can be exploited for separation; e.g. electric fields have 
been used to align nanotubes during PECVD synthesis 
1, 2] and to separate different tubes in solutions 0]. A 
detailed physical understanding of dielectric response is 
also needed to characterize optical excitations, screening 
at contacts, plasmons in nanotube arrays, and the de- 
gree of control achievable on endohedral fillings. While 
in recent years the response of SWNTs has been studied 
with tight-binding |J,|5j,|6j and first-principles approaches 
0, , MWNTs - a more common product of synthesis - 
have received much less attention due to their complexity. 
We present here a comprehensive and detailed picture 
of dielectric screening in SWNTs, MWNTs and bundles, 
using a combination of first-principles techniques and in- 
troducing an accurate classical electrostatic model that 
captures the unusual response of these materials. 

All calculations are performed using Quantum- 
ESPRESSO ;9| with the PBE approximation and ultra- 
soft pseudo-potentials in a plane-wave basis. A tetrago- 
nal unit cell is set up with periodic-boundary conditions 
in all three dimensions. A fc-point sampling grid of at 
least 30x1x1 points is used; this is sufficient to converge 
polarizabilities to within 3 significant digits. Atomic con- 
figurations are generated using an interatomic distance of 
1.42A, obtained from careful relaxation studies [lo|. Lon- 
gitudinal and transverse polarizabilities are calculated 
using density-functional perturbation theory (DFPT) 
[llj and finite-field or electric-enthalpy approaches [l2j |. 



also implemented in our Quantum-ESPRESSO code. 
Since we use periodic-boundary conditions, we effectively 
simulate a three-dimensional bulk material consisting of 
a square array of infinite parallel nanotubes. The lon- 
gitudinal dielectric response of an isolated nanotube is 
characterized by polarizability per unit length an, which 
is related to the separation-dependent bulk dielectric con- 
stant en using the relation 
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where f2 = L 2 is the cross-sectional area of the unit cell. 
From linear-response theory 0] we expect the static di- 
electric constant to depend on the gap as e(q) ~ 1 + 
(hui p /Ag) 2 which suggests via that a\\ ~ 1/A 2 ,. Our 
calculations confirm this behavior in zigzag nanotubes, 
as shown in Fig. ^ As expected, (9,0), (12,0) and (15,0) 
nanotubes have the smallest gaps and the largest a\\; the 
inverse-square dependence on the gap roughly holds over 
two orders of magnitude. Only the narrowest nanotubes 
(7,0) and (8,0) deviate from this trend. The agreement 
is particularly accurate for large-gap zigzag nanotubes 
(3n+l,0) and (3n+2,0) with n > 2. We note in pass- 
ing that for these SWNTs our first-principles results can 
be fitted well with these relations: A g as 3.3/i?o + 0.06 



and 
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20.5, with A g in eV and Rq in A. 
Previous tight-binding studies Q reported values of an 
comparable to ours, and noted a relation a\\ ~ Ro/A 2 
which we also observe for large-gap nanotubes (see Fig. 
0. For infinitely- long armchair SWNTs the longitudinal 
polarizability per unit length an diverges since there is 
no gap in the band dispersions. To get a sense of scaling 
we can approximate such nanotubes as metallic ellipsoids 
of length I and transverse radius R (I ^> R); the classical 
result is a„ ~ l 2 /[24(ln(l/R) - 1)]. For MWNTs, the lon- 
gitudinal picture remains simple: depolarization effects 
along the axis are negligible, and constituent tubes have 
very weak dielectric interactions. The total polarizability 
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TABLE I: Radius, band gap, long: 
polarizabilities (per unit length) of 
function of the chiral vector (n, m). 



itudinal and transverse 
carbon nanotubes as a 



TABLE II: Transverse polarizabilities of MWNTs. 



(n, vri) 


Ro (A)" 


A„ (eV) 


a± (A 2 ) 


an (A 2 ) 
II v / 


(7 0) 


2.74 


0.48 


6.47 


83.0 


(8 0) 


3.15 


0.57 


7.80 


104 


(9 0) 


3.58 


0.17 


9.32 


1460 


do oi 


3.95 


0.91 


10.9 


142 


(11 0) 


4.34 


0.77 


12.7 


186 


(12 0) 


4.73 


0.087 


14.3 


6140 


(13,0) 


5.09 


0.72 


16.3 


224 


(14,0) 


5.48 


0.63 


18.4 


279 


(15,0) 


5.88 


0.041 


20.3 


11100 


(16,0) 


6.27 


0.61 


22.9 


326 


(17 0) 


6.66 


0.53 


25.2 


395 


(8,0)+(17,0) 






25.8 


499 


(8,0)+(16,0) 






23.6 


427 


(4 4) 


2.71 


(0) 


6.41 


(oo) 


(5,5) 


3.40 




8.71 




(6,6) 


4.10 




11.6 




(7,7) 


4.76 




14.7 




(8,8) 


5.45 




18.1 




(9,9) 


6.12 




21.8 




(10,10) 


6.78 




26.1 




(12,12) 


8.14 




35.8 




(14,14) 


9.50 




47.2 





a Ro is the radius of the carbon backbone 



ajj ot should then simply be the sum of the polarizabili- 
ties of constituent SWNTs; this conclusion is confirmed 
by our results in Table |U 
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FIG. 1: Log scale plot of ay of zigzag nanotubes as a function 
of band gap. The dashed line has slope -2. The inset shows 
the values for large-gap SWNTs as a function of 7?o/A 2 . 



We address the characterization of the transverse di- 
electric response in two different ways. First, we calcu- 
late with DFPT the dielectric constant e±, from which 
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FIG. 2: Convergence of a± and a b j_ with respect to L for a 
(5,5) SWNT. The point at L=10.6A corresponds to a typical 
tube-tube separation in a bundle. 



the transverse polarizability a± is extracted. To study 
non-linear finite-field effects, we also obtain a± by ap- 
plying an electric field E out via a sawtooth potential, 
and computing the total induced dipole moment per unit 
length p±. In the linear regime the two approaches are 
equivalent, and we find an agreement between the two 
methods within 1%. Both DFT-based calculations take 
into account the local-field effects, that would be absent 
in tight-binding calculations. Again, since calculations 
provide us with the transverse response of a periodically- 
repeated array of nanotubes, it is necessary to remove the 
depolarization fields stemming from the periodic images. 
In principle one could use 

" V-D = ^ (2) 
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for the first and second methods respectively, while tak- 
ing the limit L — > oo for which the depolarization fields 
vanish. In practice, these persist for very large inter-tube 
separations due to the long range of electrostatic inter- 
actions between image tubes. Computation time grows 
as L 2 at a fixed energy cutoff, quickly becoming unman- 
ageable without even reaching a converged result; Fig. [21 
illustrates the slowness of this convergence. It is clear, 
however, that at large separations only electrostatic ef- 
fects are important, so we can solve this problem using 
a classical 2D Clausius-Mossotti correction [l4( relating 
the single-tube polarizability a± to the periodic bulk L- 
dependent value a b ± . The relevant conversions are 
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The values of ai obtained from © are listed in Table |U 
and plotted in Fig. |3]as a function of the square of the ef- 
fective outer radius R = Rq + 1.3A (see later discussion). 
Remarkably, transverse polarizabilitics of both metallic 
and semiconducting SWNTs lie on the same curve, which 
can be fitted by a line a± — cR 2 with slope c = 0.40. 
Thus chirality and longitudinal band structure have a 
negligible effect on the transverse dielectric response; this 
was observed in earlier calculations 0, 0, 0] and justi- 
fied with symmetry arguments in the single-particle a 
proximation Q. Recent tight-binding calculations 0, 
predict a small and systematic difference between polar- 
izabilities of metallic and semiconducting SWNTs; how- 
ever, we do not detect these differences in our DFT cal- 
culations. 
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FIG. 3: Transverse polarizabilities a± of armchair and zigzag 
nanotubes as a function of R 2 . The dashed line is the best-fit 
result of our semi-metallic shell model; the solid line a± — 
jR? corresponds to an ideal metallic cylinder. 



Periodic boundary conditions allow us to easily exam- 
ine the bulk dielectric response of nanotube bundles. We 
compute e_L and e|| of triangular and square arrays with 
inter-tube separation of ri=3. The values of en 

accurately match those computed from an of isolated 
nanotubes using (Q, thus reflecting the additive prop- 
erty of the longitudinal response. In contrast, transverse 
response of bundles depends strongly on d. Fig. [21 illus- 
trates both the benefits of using the Clausius-Mossotti 
relation (J3J) for large d, and the limitation of its applica- 
bility when d is small. Whereas the longitudinal response 
remains simple, the transverse dielectric tensor at small 
d may have sizeable anisotropic and off-diagonal contri- 
butions depending on the combined point-group symme- 
try of the nanotube and the lattice. These contributions 
vanish quickly with d and do not affect our isolated tube 
calculations. 

By applying a finite transverse field E ut we can also 
study screening inside a nanotube; we find the inner 
field Ei n to be very uniform, as shown in Fig. 0] An- 



other remarkable feature is that the screening factor 
Eout/Ei n w 4.4 ±0.1 turns out to be independent of ra- 
dius and chirality for all SWNTs. To make physical sense 
of these general results we look for a simple electrostatic 
model that would capture these traits. A solid dielec- 
tric cylinder of radius R and bulk dielectric constant e 
would have polarizability a± = ^ fqjyi? 2 , a uniform inner 
field and a screening factor E out /Ei n = (e + l)/2 inde- 
pendent of radius. This picture, however, does not corre- 
spond to a nanotube, where screening is accomplished by 
a thin layer of delocalized 7r-electrons. One could then 
treat a nanotube as a dielectric cylindrical shell of finite 
thickness. In this case the inner field remains uniform, 
but the screening factor decreases with increasing radius. 
To identify an appropriate model that incorporates all 
the observed features, we note that in general a radius- 
independent uniform inner field is produced by the sur- 
face charge density ct(</>) = do cos(^), where is the angle 
measured from the direction of E out . The dipole moment 
per unit length in this case is p± = ttctqR 2 = a±E out and 
the polarizability is a±_ = j^ t R 2 = cR 2 with c < 1/2. 
In a metallic cylinder c is 1/2 and the outer field is 
completely screened. A best fit of our ab-initio data 
for SWNTs to this model (see Fig. 0J yields the slope 
c = 0.40 and effective radius R = Rq + 1.3A larger than 
the radius of the carbon backbone i?rj, consistent with 
the finite thickness of the electronic charge density dis- 
tribution. Elementary electrostatic considerations yield a 
screening factor E out /Ei„ = r^B — 5 in good agreement 
with our finite-field calculations and previous estimates 
0, . It should be stressed that the screening properties 
of nanotubes, reflected in this model, are neither metallic 
nor insulating. This peculiarity is physically grounded in 
the fact that in a single sheet of graphite the screening 
of Coulomb interactions is anomalous due to the vanish- 
ing density of states at the Fermi points For car- 
bon nanotubes (as opposed to boron-nitride nanotubes), 
the semi-metallic nature of 7r-electrons implies that the 
screening factor is radius-invariant. 

The generalization of this model to the multi-wall case 
needs to take into account screening and electrostatic in- 
teractions between layers. Our strategy is to first solve 
exactly the general problem of N concentric dielectric 
cylindrical shells in a uniform field. We then recover pre- 
cisely the above single-layer model by treating a SWNT 
as a shell of radius R, dielectric constant e and vanish- 
ing thickness 5, and constraining these parameters by the 
scale-invariance condition 
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that guarantees that the screening factor remains inde- 
pendent of R. Modelling a general MWNT amounts to 
solving a linear system of 2N x 2N boundary-condition 
equations |l7| containing the best-fit parameters c and 
R (carried over from the single- wall case) and subject 
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in a uniform electric field. The torque on a nanotube of 
length I at an angle 9 to the field E is 



FIG. 4: Electrostatic potential for a (10,10) SWNT in an 
applied homogeneous transverse field E ou t. The electric field 
through the center slice is shown in the inset. 



to constraint (J3J. For the double- and triple- wall cases 
we find excellent agreement between this model and our 
ab-initio results (see Table |nj. We conclude that the 
present semi-metallic shell model captures all character- 
istics of the transverse dielectric response: uniform inner 
field, radius- independent screening factor in SWNTs, and 
correct a± for MWNTs. We note also that the largest 
contributions to transverse polarizabilities come from the 
outer few layers, and inner layers play a negligible role 
due to a combination of screening and their smaller radii. 
Diameter control alone thus becomes the key growth- 
parameter determining transverse response. 

The finite-field approach is also used to determine the 
range of fields for which the transverse dielectric response 
is linear. The (5, 5) nanotube exhibits precisely linear re- 
sponse with the same polarizability coefficient to within 
3 significant digits for field magnitudes of 0.05, 0.5, 5 
V/nm, the last one being greater than the experimen- 
tally attainable value. This implies that our electrostatic 
shell model of transverse response remains valid in the 
regime of large applied fields. To study the linearity of 
longitudinal response, we minimize directly the electric- 
enthalpy functional [l2| to introduce a finite longitudinal 
field while preserving periodic-boundary conditions. We 
find that the longitudinal response of the (8,0) nanotube 
becomes nonlinear by only 5% at -Em=0.5 V/nm. Nonlin- 
earity is in fact suppressed because zigzag and armchair 
(non-chiral) nanotubes are center-symmetric, so the first 
hyper-polarizability (3 vanishes by symmetry Q • To esti- 
mate the second hyper-polarizability 711 we compute po- 
larizations at several values of the field, and fit the result 
to the expression P = ot\\E + J\\E 3 . We obtain aii = 106 
A 2 (in agreement with the DFPT result in Table and 
7|l = 3.1 x 10 7 in atomic units. 

We now turn to the question of alignment of nanotubes 



t = |p x E| = 1 1 
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The longitudinal and transverse polarizabilities compete 
with each other, but our results imply that an > aj_ 
in all nanotubes, much more so in metallic and small- 
gap semiconducting nanotubes. Indeed, for all nanotubes 
a± < \R 2 whereas for large-gap SWNTs ay > 8.2i? 2 , 
and for MWNTs an is additive while a± is not. So nan- 
otubes of all types will align with the electric field, but 
by tuning the value of the field during PECVD growth 
it may be possible to selectively grow highly polarizable 
(e.g. metallic) tubes. 

There have also been attempts to separate semicon- 
ducting and metallic nanotubes in solution using inhomo- 
geneous electric fields Q. A polarized nanotube aligned 
with the field will be pulled in the direction of or against 
the field gradient, depending on its effective dielectric 
constant en relative to that of the solvent e s . Assuming 
no solvent inside the nanotube, and approximating it by a 
solid dielectric cylinder of radius Rq, we obtain from our 
values of a\\ an effective e|| = 1 + Aa^/R^ k, 30 for large- 
gap semiconducting SWNTs and obviously much larger 
values for small-gap and metallic tubes. This result is 
consistent with findings that only metallic SWNTs are 
observed deposited on the electrodes in water (e s = 80), 
whereas all nanotubes are drawn towards the electrodes 
in isopropyl alcohol (e s = 18). 

In summary, we studied in detail the dielectric prop- 
erties of isolated and bundled SWNTs and MWNTs. In 
SWNTs, the longitudinal response is controlled by the 
band gap, while the transverse response is sensitive only 
to the effective radius. In bundles and MWNTs longitu- 
dinal response is additive, while the transverse response 
in MWNTs is dominated by the outer few layers. We pre- 
sented an accurate scale-invariant electrostatic model of 
transverse response, which is intermediate between that 
of a metal and an insulator. The authors would like 
to thank L. S. Levitov for valuable suggestions. This 
work was supported by NSF-NIRT DMR-0304019 and 
the Singapore-MIT alliance. 
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